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ABSTRACT 

Context. Previous research has indicated that high amounts of ozone (O3) and oxygen (O2) may be pro- 
duced abiotically in atmospheres with high concentrations of CO2. The abiotic production of these two 
gases, which are also characteristic of photosynthetic life processes, could pose a potential "false-positive" 
for remote-sensing detection of life on planets around other stars. We show here that such false positives are 
unlikely on any planet that possesses abundant liquid water, as rainout of oxidized species onto a reduced 
planetary surface should ensure that atmospheric H2 concentrations remain relatively high, and that O2 and 
O3 remain low. 

Aims. To determine the amount of O3 and O2 formed in a high CO2 atmosphere for a habitable planet with- 
out life. 

Methods. We use a photochemical model that considers hydrogen (H2) escape and a detailed hydrogen 
balance to calculate the O2 and O3 formed on planets with 0.2 of CO2 around the Sun, and 0.02, 0.2 and 
2 bars of CO2 around a young Sun-like star with higher UV radiation. The concentrations obtained by the 
photochemical model were used as input in a radiative transfer model that calculated the spectra of the 
modeled planets. 

Results. The O3 and O2 concentrations in the simulated planets are extremely small, and unlikely to pro- 
duce a detectable signature in the spectra of those planets. 

Conclusions. With a balanced hydrogen budget, and for planets with an active hydrological cycle, abi- 
otic formation of O2 and O3 is unlikely to create a possible false positive for life detection in either the 
visible/near-infrared or mid-infrared wavelength regimes. 

Key words, planetary atmospheres - biosignatures - Terrestrial Planet Finder 



1. Introduction 

Most of the O2 in Earth's present atmosphere is thought to have been produced by oxygenic 
photosynthesis, followed by burial of organic carbon in marine sediments (Cloud [1972] Walker 
IT9771 Hollan d 1978lfT984l[2002b . Prior to the origin of life, and of Oa-producing life in particular, 
atmospheric O2 mixing ratios are thought to have been very low, 10 by volume, or ~ 10"'^ 
PAL at the surface (Walker [T9771 Kasting et aL [T979l Kasting [HSU Kasting & Cading |2003]l 
Here, 'PAL' means 'times the Present Atmospheric Level', which is 21 percent by volume, or 
0.21 bars. O2 concentrations in the upper atmosphere of the early Earth could have been much 
higher, up to ~ 10 by volume (Kasting and Catling l2003l l. as a consequence of photolysis 
of CO2, followed by recombination of O atoms to make O2. A small amount of ozone, O3, 
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could have been formed from this O2, but not enough to provide an effective shield against 
solar UV radiation (Levine et al. 119791 Kasting and Donahue 119801 Kasting 119931 and refs. 
therein; Cockell and Horneck. 1200 11 1. Early life therefore probably evolved in an anoxic, high- 
UV environment. 

While the question of abiotic O2 and O3 levels has historically been of interest to geologists 
and biologists, it has recently become an important issue for astronomers as well. Within the next 
15-20 years, space-based telescopes, such as the two Terrestrial Planet Finder (TPF) missions 
planned by NASA (http://planetquest.jpl.nasa.gov/TPF/tpf_c_team.cfm) and ESA's Darwin mis- 
sion |Ettp://ast.star.r l. ac.uk/darwin/), will hopefully search for Earth-sized planets around other 
stars and attempt to obtain spectra of their atmospheres and surfaces. TPF-C (a visible/near-IR 
coronagraph) will be sensitive to the 0.76-/im band of O2, while TPF-I and Darwin (both thermal- 
IR interferometers) will be sensitive to the 9.6-/vmband of O3. Both of these gases are considered 
to be possible biomarker compounds (Owen 1980; Angel , 1986: Leger et al. 119931 ). The question 
as to whether they can be produced abiotically is therefore of great potential relevance to the 
interpretation of data from these future missions. 

Situations in which either O2 or O3, or both, might accumulate abiotically have been iden- 
tified by a number of different authors. Some of these situations are legitimate "false positives" 
for Ufe, while others might be misleading. Two situations, in particular, appear capable of pro- 
ducing "false positive" signals. The first is a runaway greenhouse planet, like early Venus, on 
which large amounts of hydrogen escape from a hot, moist atmosphere (Kasting [l997i Schindler 
& Kasting |2000l ). Because the hydrogen originates from H2O, oxygen is left behind. The escape 
of a terrestrial ocean equivalent of hydrogen, unaccompanied by oxygen sinks, could leave an 
atmosphere containing ~ 240 bars of O2 (Kasting 1997 ). A second conceivable "false positive" 
is a frozen. Mars-like planet that is large enough to retain heavy gases but too small to maintain 
volcanic outgassing (Kasting 1997; Schindler & Kasting 2000). The frozen surface would in- 
hibit the loss of oxygen by reaction with reduced minerals, whereas the lack of outgassing would 
eliminate reaction with reduced volcanic gases (primarily H2) as an oxygen sink. The martian 
atmosphere contains 0. 1 % O2 and would likely have even more if the planet were slightly larger 
so that it did not lose oxygen to space by nonthermal loss mechanisms (McElroy and Donahue 
[19721) . 

Both of the "false positives" mentioned above apply to planets that lie outside of the liquid 
water habitable zone around their parent star The boundaries of this zone can be estimated to 
first order from climate models (Kasting et al. 119931 ), and it should be possible to determine 
observationally whether the planet's distance from its parent star falls outside these limits. Such a 
determination would not be definitive, as the theoretical limits of the habitable zone are uncertain, 
mostly because of the difficulty in simulating clouds. (The Kasting et al. 1 19931 model effectively 
puts the cloud layer at the planet's surface, thereby ignoring all cloud feedbacks.) Planets outside 
the habitable zone may also be distinguishable by the absence, or near absence, of gaseous H2O 
in their spectra, although water-rich planets near the inner edge might not obey this rule, and one 
should remain suspicious of O2 signals on such bodies. Still, the most important false positive 
issue is whether or not planets within the habitable zone could build up O2 or O3 abiotically. 

Calculations predicting high, or relatively high, abiotic O2 and O3 concentrations have ap- 
peared sporadically in the literature over the past 40 years. All authors have realized that photol- 
ysis of H2O, followed by escape of hydrogen to space, is a net source of oxygen. Berkner and 
Marshall ( ll964lll965|[li66l[T^67b estimated abiotic O2 concentrations of 10"* - 10"^ PAL based 
on how much O2 buildup was needed to block out the UV radiation that dissociates H2O. This 
is probably too little O2 to be detectable spectroscopically by a telescope like TPF-C; however, 
the associated ozone layer could conceivably be detected by a telescope like TPF-I or Darwin 
(Segura et al. 2003). In a study a few years later, Brinkman (1969) found abiotic O2 concen- 
trations of up to 0.27 PAL. If correct, such an atmosphere would produce an absorption feature 
at 0.76 yum about half as strong as that of Earth (Des Marais, et al. 120021) . Brinkmann obtained 
this high value because he assumed that precisely one-tenth of the H atoms produced by H2O 
photolysis escaped - a fraction that we now know is much too high. 

All of these early studies were performed before the factors controlling hydrogen escape 
from Earth's atmosphere were well understood. Hunten (1973 ) showed that the H escape rate 
from Earth (and from Saturn's moon. Titan) is limited by the rate at which hydrogen can diffuse 
upwards through the homopause. (The homopause, near 90 km in Earth's atmosphere, is the al- 
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titude above which light gases begin to separate out from heavy gases. Equivalently, it marks the 
transition from turbulent vertical transport to molecular diffusion.) The resulting escape rate is 
termed the "diffusion-limited flux". Hydrogen can escape more slowly than the diffusion limit 
because of energy limitations higher up (see, e.g., Tian et al. 2005), but it cannot escape more 
rapidly. Walker (|1977| l applied this concept to the early Earth and showed that it implies ex- 
tremely low ground-level prebiotic O2 concentrations, of the order of 10 PAL. This same 
reasoning applies to vertically resolved atmospheric photochemical models (Kasting et al. 1 19791 
Kasting [l993l Kasting and Catling l2003l l. although the O2 concentration in the stratosphere can 
be much higher (10^^), as mentioned earlier. High stratospheric O2 concentrations are favored 
by high CO2 levels and high solar UV fluxes (Canuto et al. 119821 11983i It should be pointed 
out that the extremly high UV fluxes studied by Canuto et al. are unlikely to apply to the early 
Earth or to any analogous planet, as the T-Tauri phase of stellar evolution was probably over long 
before the Earth was fully formed. 

Although most recent models predict low abiotic O2 concentrations for planets within the 
habitable zone - contradictory results also appear. This happens, for example, in some calcu- 
lations performed by Selsis et al. (2002 ) for a variety of Earth-like and Mars-like planets. The 
Earth-like planets all had significant volcanic sinks for oxygen, and so none of these cases pro- 
duced high O2 or O3. However, their "early Mars-type" planet (Case B2), which assumed zero 
volcanic outgassing, exhibited 0. 1 PAL of O2 and a "super" ozone layer with a column depth of 
~ 0.7 atm cm - roughly twice that of modern Earth. This model simulated a warm, humid, 1-bar 
CO2 atmosphere underlain by an ocean. If this calculation were valid, then such a planet would 
represent another possible "false positive" for life. Selsis et al. argued that the O3 9.6-fim band 
would be obscured by the neighboring 9.4- and 10.4-jum hot bands of CO2, and thus would not 
constitute a false positive for Darwin (or TPF-I). However, the 0.76-;um band of O2 would be 
quite prominent in this atmosphere, and so it would remain an issue for TPF-C. 

We argue here that the Selsis et al. calculation is incorrect, or at least over-simplified, and that 
this case is not a likely "false positive". More precisely, we suggest that abiotic O2 buildup is only 
likely to occur on dry or frozen planets, and such planets could be differentiated from warmer 
Earth-like planets by examining their spectra. The reason that O2 reached such high levels in 
the Selsis et al. calculation is that the authors failed to consider the effect of rainout of oxidized 
(or reduced) species on the atmospheric hydrogen budget. Photochemical production of oxidized 
species such as hydrogen peroxide, H2O2, followed by their reaction with reduced species in 
the crust or in seawater, should generate a net source of H2 (Kasting et al. '1984a; Kasting and 
Brown [19981 Kasting and Catling 2003J . This H2 source (and O2 sink) was left out of the Selsis 
et al. model, allowing for a much higher abiotic O2 concentration. The next section describes the 
atmospheric hydrogen budget, along with the rest of our photochemical model, in more detail. We 
then use our model to simulate various Earth-like planets, with and without volcanic outgassing, 
and under high and low stellar UV flux, and we calculate abiotic O2 and O3 concentrations and 
accompanying visible and thermal-IR spectra. Our goal is to determine whether any additional 
"false positives" for O2 producing life could exist in high-C02 atmospheres. 

2. Model description 

2.1. The photochemical model 

The 1-D photochemical model used in our study was developed for high-C02/high-CH4 ter- 
restrial atmospheres by Pavlov et al. (120011 1 and was subsequently modified by Kharecha et al. 
(I2OO5I 1. The model simulates an anoxic atmosphere composed of 0.8 bar of N2 and variable 
amounts of CO2. We did not simulate a pure-C02, "Mars-like" atmosphere explicitly; however, 
such an atmosphere is similar in principle to our high-C02 cases. Our photochemical model con- 
tained 73 chemical species involved in 359 reactions and spanned the region from the planetary 
surface up to 64 km in 1-km steps. The solar zenith angle was fixed at 50°, and a two-stream 
approach was used for the radiative transfer. The continuity equation was solved at each height 
for each of the long-lived species, including transport by eddy and molecular diffusion. The 
combined equations were cast in centered finite difference form. Boundary conditions for each 
species were applied at the top and botton of the model atmosphere, and the resulting set of cou- 
pled differential equations was integrated to steady state using the reverse Euler method. Two 
different stellar UV fluxes were used, as described below. The eddy diffusion profile for most 
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cases was assumed to be the same as the one measured on present Earth (Massie & Hunten 
119811 1. with modifications for atmospheres in which the total surface pressure exceeded one bar. 
All of the simulated planets were assumed to be devoid of life; hence, none of the compounds in 
the atmosphere was considered to have a biological source. 

As a reference, or "standard", case we assumed a 1-bar, cloudless, "early-Earth" type atmo- 
sphere with 0.2 bar of CO2. This is roughly the amount of CO2 needed to compensate for a 30% 
reduction in solar luminosity, appropriate for a time near 4.5 Gyr ago. (KastingTWT). We did 
not perform coupled climate model calculations for this case; rather, we simply assumed that the 
temperature decreased with height from 278 K at the surface to 180 K at 12.5 km, following a 
moist adiabat. Above that height, the atmosphere was assumed to be isothermal (Fig.[T]). This is 
roughly consistent with the predictions of climate models (Kasting 1990). It is difficult for us to 
improve significantly on this assumption, as our own climate model (Pavlov et al. 120001 1 is not 
particularly accurate up in the Doppler-broadened stratosphere and mesosphere. The assumed 
solar UV flux in this model was the same as for Earth today. We also performed simulations for 
higher UV fluxes. As discussed further below, stellar UV fluxes at wavelengths less than ~200 nm 
originate from a star's chromosphere and are, if anything, anti-correlated with the flux at longer 
wavelengths. Stars tend to brighten in the visible as they age; however, their short- wavelength 
UV emissions decrease with age in most stars because their spin rate slows and they exhibit less 
chromospheric activity. Longer wavelengths are included in the model, but they mostly affect the 
temperature profile, which has been fixed for the planets simulated here. 

We also did calculations for early Earth-type planets with CO2 partial pressures of 0.02, 0.2, 
and 2 bars and higher UV radiation than the present solar flux. For the 0.02 -bar case, we kept the 
surface temperature and pressure the same as for the standard model. (This implies slightly more 
N2 and either a higher incident solar flux or additional unspecified greenhouse warming.) For the 
2-bar CO2 case, the total surface pressure was 2.9 bars (- pC02 + pN2 + PH2O). The surface 
temperature, 317 K, and the vertical temperature profile for this atmosphere were calculated by 
a climate model (Kasting |l9901 l. In this latter case, the eddy diffusion profile was readjusted to 
account for the increased troposphere thickness. Following Kasting (1990") we shifted the eddy 
diffusion profile upward by 4 km for the 2 bar case, and assigned an eddy diffusion coeflicient 
of 10^ cm^ s ' to the 4 lowermost grid points (Fig|2]). Although the cases that we studied are 
not entirely self-consistent, we chose them because they span a plausible range of atmospheric 
composition for abiotic, early Earth-type planets (Kasting 1990). 

Volcanic outgassing of H2 and other reduced species was included in most, but not all, of the 
model calculations. For H2, we used the present volcanic outgassing rate ~ 5 x 10'^ mol yr"\ or 
~ 2x10'" cm"^ s ' (Holland l2002l ). Sleep and Bird (2007) argue for an H2 outgassing rate that is 
lower by a factor of 10. Our calculations span this range because we consider zero outgassing as 
well. The methane flux for the standard case was set to the estimated modern non-biogenic flux 
on hydrothermal vents. Formerly, the biotic-to-abiotic CH4 ratio was estimated to be ~ 300 : 1 
(Kasting & Catfing |20031 l. New measurements of methane dissolved in the fluids of the Lost City 
vent field (Kellev l20051 l indicate that the abiotic methane flux is higher than previously thought 
by a factor of 10; therefore, the biotic-to-abiotic ratio may be ~ 30 : 1. The present biological 
methane flux is 535 Tg CH4/year or 3.3 x 10'^ mol/yr (Houghton et al. 119951 1. so the assumed 
abiotic methane flux was 1 x 10'2 mol/yr, or 3.7 x 10*^ cm -s . The actual abiotic methane flux 
may be lower than this value if the CH4 in the vent fluids is produced by methanogens living deep 
within the vent systems. An (implicit) NH3 flux of 3.87 x 10"^ cm"^ s"' was included because 
the standard model contained 10 ppbv of NH3. This was simply for convenience, because NH3 
was included in the species list and, hence, its concentration must be non-zero. For the "no- 
outgassing" case, the methane and H2 fluxes were set equal to zero, but volcanic outgassing of 
SO2 and a tiny amount of NH3 was retained in order to avoid numerical problems associated with 
disappearing species within the code. The implications of this assumption are discussed below. 

Our model also includes the volcanic outgassing of SO2 at its present estimated release rate 
3.5 X 10^ cm ' s~' (Kasting 19901 ) along with formation of NO, O2 and CO by lightning in the 
troposphere, as described by Kasting (ll9901 l.Rainout rates for soluble species are calculated using 
the method of Giorgi & Chameides (119851 1. These rainout terms are added to the normal chem- 
ical loss rate for each species. Sulfur gases are removed from the model atmosphere by rainout 
and surface deposition and by conversion into particulate sulfur and elemental sulfur, followed 
by rainout and surface deposition of the particles. Rainout lifetimes of high soluble gases (e.g. 
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Fig. 1. a) Temperature profiles for the I0W-CO2 atmospheres (0.2 and 0.02 bars of CO2) and the 
high CO2 planet (2 bars of CO2. b) Water profiles for the simulated planetary atmospheres around 
EK Dra. 



H2SO4) and of particles was approximately five days. Less soluble gases (SO2, H2S) have longer 
lifetimes against rainout, according to this model. The lifetimes assumed here are appropiate for 
the modern atmosphere, we assume that the hydrological cycle and, hence, rainout lifetime was 
similar on early Earth. This assumption seems reasonable, since evaporation (and, thus, precipi- 
tation) rates are ultimately controlled by solar heating rates, not by temperature (Holland. ll978l l. 
Furthermore, the rainout rate of highly soluble gases is determinated by how often it rains, no 
how much it rains. Soluble species are also removed by direct deposition at the lower boundary, 
simulating uptake by the ocean. The downward flux of a chemical species i is equal to its number 
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Fig. 2. Eddy diffusion coefficient (^edd) profiles used in our simulations. Present Earth fcgdd is the 
one measured by Massie & Hunten (I1981I I. the other two profiles are modified versions of this 
coefficient for the 2 bar of CO2 atmosphere and more pausible kedd for early Earth (see Results). 

density, n,, multiplied by a deposition velocity, Vdep. Deposition velocities for different gases are 
taken from Slinn et al. (,1978j and Lee & Schwartz (1981) . Values range from ~ 0.02 cm/s for 
less reactive species (e.g., H2O2) to 1 cm/s for the most reactive species (e.g., OH). The upper 
limit on Vdep is set by diffusion through the turbulent atmospheric boundary layer, assuming that 
molecules are absorbed by the surface each time they collide with it. 

One may ask how our photochemical model achieves steady state when species are continu- 
ally injected into the atmosphere from volcanoes. For hydrogen and sulfur, the answer is clear: 
they are removed either by escape to space (for H) or by rainout from the atmosphere. For CH4 
and NH3, the answer is less obvious. These species are oxidized photochemically to CO2 and 
N2, respectively. As the CO2 and N2 concentrations are held constant in the model atmosphere, 
these species are implicitly lost by reactions of CO2 with surface minerals (to form carbonates) 
and by unspecified mechanisms of N2 loss (most plausibly biological nitrogen fixation). 

2.2. The atmospheric hydrogen budget 

A key feature of our model that makes it useful for this study is its ability to keep track of 
the atmospheric hydrogen budget, or redox budget. The basic principle behind this budgeting 
scheme is simple: when one species is oxidized, another species must be reduced, and vice versa. 
Alternatively, one could rephrase this principle as requiring that the electron budget of the model 
atmosphere be balanced. 

In practice, it is easiest to keep track of redox balance in terms of the abundance of molecular 
hydrogen, H2. Following Kasting and Brown d 19981 1 we define "redox-neutral" species: H2O (for 
H), N2 (for N), CO2 (for C), and SO2 (for S). All other species are assigned redox coefficients 
relative to these gases by determining how much H2 is produced or consumed during their forma- 
tion from redox neutral species. For example, formation of hydrogen peroxide can be expressed 
as: 

2H2O H2O2 + H2 

The coefficient of H2O2 in the hydrogen budget is therefore - 1 , meaning that when one H2O2 
molecule is rained out, one H2 molecule is produced. Similarly, formation of elemental sulfur, 
Sg, can be written as: 
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8SO2 + I6H2 ^ Sg + I6H2O 

Hence, rainout of Sg particles consumes 16 H2 molecules (per Sg molecule), so the redox 
coefficient of Sg is +16. Volcanic outgassing of species other than H2 itself is treated in the 
same manner, i.e., it adds to the total hydrogen outgassing rate according to the appropriate 
stoichiometry. As mentioned earlier, the standard model includes modest outgassing of CH4 and 
NH3, in addition to H2. Thus, the total hydrogen outgassing rate can be written as: 
<l>voic(H2) = <1)(H2) + 1.5<D(NH3) + 4(D(CH4) 

where the terms on the right represent the fluxes of the individual reduced species multiplied by 
their stoichiometric coefficients in the redox budget 

With these definitions in place, the hydrogen budget can be written as: 

<I>volc(H2) + <l>rain(Ox) = (De.sc(H2) + <l)rain(Red) (1) 

Here, the terms <l>,ain(Ox) and <l>rain(Red) represent the net flux of photochemically produced 
oxidants and reductants, respectively, from the atmosphere to the ocean. This flux includes both 
rainout and surface deposition. Eq. ^ is diagnostic, not prognostic; hence, it provides a good 
check both on redox balance and on the photochemical scheme. (This equation will not balance 
if any of the chemical reactions are not balanced.) In our model, Eq. ([T]l balances to within one 
part in 10^ for the standard case and one part in 10^ for the other simulated planets (Table [1]). 

As discussed earlier, we assume that the hydrogen escape is limited by the H2 diffusion rate 
through the homopause. The diffusion-limited escape flux is given by Hunten ( 1973i ) and Walker 
( fT977b : 

(2) 

s 2.5 X 10 /tot (molecules cm ^ s ') 

where b is an average binary diffusion coefficient for the diffusion of H and H2 in nitrogen, 
H is the atmospheric scale height (RT/g), and /tot is the sum of the mixing ratios of all the 
hydrogen-containing species, weighted by the number of hydrogen atoms they contain. As we 
have expressed the hydrogen budget in terms of H2 molecules, we do the same for /ot: 

/tot = /(H2) + 0.5/(H) + /(H2O) + . . . (3) 

Because the stratosphere is cold and dry in all of our simulations, the contribution of H2O to 
ftot(H2) is negligible. When this is not the case (as in the runaway greenhouse example mentioned 
earlier), high abiotic O2 concentrations are possible. The actual hydrogen escape rate may be 
lower than predicted by Eq. (|2]l However, that is not a concern for this paper, as slower H escape 
will lead to higher atmospheric H2 concentrations, and hence lower abiotic O2 and O3. When 
trying to identify "false positives" for life, the conservative approach is to assume the maximum, 
diffusion-limited, hydrogen escape rate. 



2.3. The input Stellar spectra 

Some of the stars observed by the TPF missions and by Darwin may be young stars with UV 
fluxes higher than that of the Sun. High stellar UV should increase the rate of CO2 photolysis, and 
thus enhance high-altitude concentrations of O2, and possibly O3, on planets orbiting such stars. 
High atmospheric CO2 concentrations on young, early Earth-type planets would presumably 
enhance production of O2 and O3. To determine whether this process might produce enough 
abiotic O2 and O3 to be detected remotely, we simulated planets with high-C02 atmospheres 
around a star with a much higher UV flux than our Sun. EK Draconis was chosen from a suite of 
solar-type stars studied by Ribas et al. (2005) because it has the highest UV flux of the six stars 
studied. EK Dra itself may be an unlikely candidate for study by TPF or Darwin because of its 
extremely young age (see Discussion). However, our goal here was to study a planet exposed to 
a maximal UV radiation flux to determine if this environment favors abiotic formation of O2 and 
O3. EK Dra (HD 129 333) is a G1.5 V spectral type member of the Pleiades moving group. It 
has an estimated age of ~ 0. 1 Gyr, making it a good proxy for the young Sun (Dorren & Guinan 
[T9941 Strassmeier & Rice 17998] Montes et al. [20011 Ribas et al. i2C)051 l. Dorren & Guinan (fr994] l 
state that it is the "only known star closely matching the zero age main sequence Sun observable 
with lUE" (the International Ultraviolet Explorer.) 
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Fig. 3. UV fluxes for the Sun and EK Dra at the top of the atmosphere of a habitable planet at 1 
AU equivalent distance. 



An absolute stellar spectral energy distribution was assembled from a Kurucz ( 119931 ) photo- 
spheric spectrum with T^ff - 5765 K, surface gravity \ogg = 4.61 g cm"^, and solar metallicity. 
(Gray et al. (1993) derived [M/H] -0.04 + 0.09 dex). The spectrum was absolutely normalized 
at 21 observed photometric points between 0.38 and 12 yum, as described by Cohen et al. (I2003I I. 
Nineteen empirical UV spectra (17 SW, 2 LW) were recalibrated (Massa & Fitzpatrick l2000l) and 
combined using the overlap. From 0.335 to 160 fim the spectrum is the normalized photosphere. 
Below 0.335 fim it is an 11 -point boxcar-smoothed version of the empirical UV data. 

The Lya flux has not been measured for EK Dra. Therefore, we estimated its flux using the 
empirical expression derived by Ribas et al. (2005 , Eq. (2)). At the top of the planet's atmosphere, 
the Lya flux was estimated to be 105.8 ergs cm"^ s ' . By comparison, the Lya flux from our Sun 
varies from ~ 6 to 10 ergs cm"^ s"' during one activity cycle at 1 AU (Emerich et al. |2005l l. 
The simulated planets around EK Dra were positioned so as to receive the same integrated solar 
flux as present Earth (1375 W m"^). This condition is satisfied if the planet's orbital radius is 
0.97 AU. (EK Dra is slightly less luminous than the present Sun.) The UV flux at the top of the 
simulated planet's atmosphere is shown in Fig.|3] Earth's UV flux is also shown for comparison. 
The solar flux used in the code comes from the World Meteorological Organization (I1985I I 



2.4. Radiative transfer model 

To generate the spectra shown in this paper, the pressure and temperature profiles and calculated 
mixing ratios from the coupled climate-chemical model are used as input to a line-by-line radia- 
tive transfer model, which generates angle-dependent synthetic radiance spectra. As the first step 
in this process the temperature and pressure profiles, and the calculated mixing ratios for H2O, 
O2, O3, CO2, and CH4 are input to a model (LBLABC: Meadows and Crisp 1996) which gen- 
erates line-by-line monochromatic gas absorption coefficients. LBLABC employs several nested 
spectral grids that completely resolve the narrow cores of individual gas absorption lines at all 
atmospheric levels, and include their contributions at very large distances from the line center 
(e.g. 500 to 1000 cm ' from the line centers.) 

This model was originally developed for simulating radiative processes in the deep atmo- 
sphere of Venus (Meadows and Crisp , 1996) and for highly precise simulation of radiances and 
heating rates in the Earth's atmosphere (Crisp 119971 ). In addition, since we are working with 
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relatively massive CO2 atmospheres, we explicitly included the effects of collisionally-induced 
absorption for CO2. Our current models include CO2 continuum absorption at all wavelengths 
where it has been measured or predicted. At wavelengths between 33 and 250 fim, we use the 
pressure induced rotational absorption (PIA) formulation described by Gruszka and Borysow 
(|1996l|1997lll998l l. At wavelengths between 4 and 10 fim, we simulate the CO2 PIA using the 
formulation described by Kasting et al. (|1984bl l. At near-IR wavelengths, we used the empirical 
pressure-induced absorption coefficients from Moore (|19711 l. 

LBLABC simulates the line shape function differently in the line-center and far-wing regions. 
For line-center distances less than 40 Doppler halfwidths, a Rautian line shape is used. This 
line shape incorporates Doppler broadening, collisional (Lorentzian) broadening, and collisional 
(Dicke) narrowing (c.f. Goody and Yung |1989] l. At greater distances, a Van Vleck-Weisskopf 
profile is used for all gases except for H2O, CO2, and CO. The super-Lorentzian behavior of the 
far wings of H2O lines, which has been attributed to the finite duration of collisions, is param- 
eterized by multiplying the Van Vleck-Weisskopf profile by a wavenumber-dependent X factor 
(Clough et al. |19891 l. This profile is adequate to account for the water vapor continuum absorption 
seen at thermal infrared wavelengths, as well as the weak continuum absorption seen throughout 
much of the visible and near-infrared spectrum (Crisp [T997] l. The sub-Lorentzian behavior of the 
CO2 line profile is primarily a consequence of a vibration-rotation energy redistribution process 
called collisional line mixing. Direct ab initio methods now exist for computing these effects 
(c.f. Hartmann and Boulet. [T991l l. but these methods do not yet provide the accuracy or numer- 
ical efficiency required for routine use. LBLABC currently employs a simple, semi-empirical 
algorithm to correct for the effects of line mixing. For strongly mixed Q-branch lines, we use 
the first order correction described by Rosenkranz ( |1988l l. With the line mixing correction and 
the use of a sub-Lorentzian Chi-factor (c.f. Fig 2 of Meadows and Crisp 1996), these methods 
provide adequate accuracy for modeling the CO2 absorption at pressures and temperatures as 
high as those seen in the deep atmosphere of Venus (90-bars), and should provide the accuracy 
needed for the calculations proposed here. 

These products, along with the atmospheric properties and constituent mixing ratios were 
used as input into our Une-by-line radiative transfer model, SMART_SPECTRA, which was de- 
veloped by the Virtual Planetary Laboratory for use with generalized terrestrial planet climate 
models, and is an updated version of the Spectral Mapping Atmospheric Radiative Transfer 
model (SMART; Meadows & Crisp [19961 Crisp [T997] l developed by D. Crisp. For this appli- 
cation, SMART_SPECTRA was used to calculate a high resolution spectrum of each planet for 
a solar zenith angle of 60°, which approximates the average illumination observed in a planetary 
disk-average. All spectra were generated over an ocean surface, which provides a neutral back- 
ground for understanding the atmospheric changes in the models. The ocean surface, however, 
is extremely dark; thus, the reflectivities shown here will typically be much lower than would 
be expected for the present Earth disk-average, which would also include the higher-reflectivity 
contributions from continents and clouds. For the visible and near-infrared spectra shown here, 
the albedo is computed as the ratio of the integrated upward flux over the hemisphere, divided by 
71 steradians, to approximate the mean radiance seen in the disk-average, and then divided by the 
solar insolation at the top of the atmosphere. For the mid-infrared spectra, the radiance shown is 
the integrated upward flux over the hemisphere, divided by tt steradians. 

3. Results 

3.1. Early Earth standard atmosphere 

The vertical profiles of O2 and O3 obtained in our "standard" atmosphere (0.2 bars CO2 planet 
around the Sun) are presented in Fig. S] and the associated columns depths are summarized in 
Tabled 

To illustrate the role of volcanic outgassing of reduced gases in determining O2 and O3 con- 
centrations, we performed a simulation with H2 and CH4 outgassing turned off. As expected, 
more O2 and O3 were formed in the no-outgassing case. The calculated increases relative to the 
standard case were: ~4 times for O2 and a factor of ~ 40 for O3 (Table |2]i. The latter increase 
is, of course, large in a relative sense; however, the absolute ozone column density for the no- 
outgassing case is still < 0.1% of the mean ozone column density on present Earth. The reason 
why the calculated increases in O2 and O3 were so modest is that H2 did not disappear from the 
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Table 1. Hydrogen budget for planets around the Sun and EK Dra (high UV) with 0.2 bars of 
CO2. 







14. 

n2 


budget contribution" 




LUClilCiCllL 


^ 1" Q Ti n r/i 
OLallQdlU 


No outgassing 


rlign U V 


Reduced species 










\j 
rl 


U.J 


H-.OO X iU 


on \y 1 n6 
S.oU X iU 


Z.4 / X lU 




1 


J.JO X lU 


3.1)2 X 11) 


Z.U4 X lU 




1 c 
V.J 


1 Al \y 1 
i.4Z X iU 


4.4j X tU 


Z.IK X lU 




9 

z 


J. 00 X lU 


X tU 


y.io X lu 


ri20 


Q 

J 


Z. 1 D X lU 


3.84 X it) 


1. /V X lU 


14C 


Z.J 


Q /in w 1 n7 
0.4U X lU 


3.y / X It) 


0.4U X lU 


14CM 


1 c 
L.J 


/ .jO X lU 


3.49 X 10 


4.JJ X lU 


M14 

iNrl2 


1 
i 


o.UZ X iU 


0. / J X lU 


y.U4 X lU 


M T4 


1 c 
L.J 


1 OA \y 1 
i.Z4 X iU 


8.10 X lU 


1 An \ r\l 
1.4/ X lU 


NT 1 1 


Z 


i.jU X iU 


3.33 X 10 


1 vy 1 n8 
l.Zj X lU 






4. 4 X lU 


1 on \y 1 n8 

j.sy X it) 


y.OJ X lU 


Oxidized species 










U 


1 

— 1 


Z.UU X lU 


3.91 X 10 


O.OJ X lU 


140 


— l.J 


A ^A V 1 

O.JO X iU 


0. /3 X 10 


0.9/ X lU 




— 1 


/.Zo X lU 


2.08 X 10 


1.1 J X lU 


14MM 


—U.J 


Z.Oj X lU 


2.96 X 10* 


3.47 X 10^ 


H2SO4 


-1 


8.64 X 10^ 


6.97 X 10^ 


2.21 X 10^ 


SO4 aerosol 


-1 


8.47 X 10' 


1 00 1 nS 

1.58 X 10 


1.33 X 10 






1 1 Q \/ 1 n8 
1 . IV X lU 


5.31 X 10* 


2.01 X 10* 


volcanism contribution 










H2 


1 


2.00 X 10'" 


0.0 


2.00 X 10'" 


NH3 


1.5 


3.87 X 10' 


5.27 X 10"^ 


4.21 X 10' 


CH4 


4 


1.80 X 10'" 


0.0 


1.80 X 10'" 


*volc(H2) 




4.19 X 10'" 


5.27 X 10"^ 


4.22 X 10'" 


Oesc(H2) 




3.76 X 10'" 


1.47 X 10* 


3.26 X 10'" 


H2 balance 




3.67 X 10' 


6.02 X 10^ 


3.53 X 10^ 



" In molecules cm - cm ' . 



Table 2. O3 and O2 column depths (cm ^) and H2 and CH4 volume mixing ratios for high CO2 
terrestrial atmospheres. Present Earth values are shown for comparison. 
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Fig. 4. 02 and H2 mixing ratios (a) and O3 number densities (b) for the standard and no outgassing 
cases. 

atmosphere even though volcanic outgassing was turned off (Fig.|4^). The H2 remaining in this 
no-outgassing case is produced by atmospheric photochemistry, followed by rainout of oxidized 
species from the model atmosphere (Table[Tll- 

A set of sensitivity tests was performed in order to quantify some of the uncertainties in our 
model. Most of these test were done for the "no outgassing" case, which maximizes the amount 
of abiotic O2 and O3. Five specific modifications were tested in the model, the results of which 
are listed in Table |2] First, we doubled the third body reaction rates for O3 and O2 formation 
to simulate the possible effects of having CO2 as a third body. This caused a modest (~ 30%) 
increase in the column depth of O3 and actually decreased the column depth of O2. Second, 
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we eliminated the water absorption for wavelengths greater than 195 nm where the absorption 
coefficients have been extrapolated from measurements at shorter wavelengths. This caused a 2- 
fold increase in both O2 and O3, presumably due to the slower recombination of O2 and H2 when 
the abundance of H and OH radicals (from H2O photolysis) is reduced. Third, the minimum of 
the eddy diffusion coefficient (feedd) was increased by 3 (Fig. |2]i in order to simulate the effect 
of no having an inversion layer due to the lack of O3 heating in the atmosphere as explained by 
Kasting et al. (1979). This decreased O2 by about 40% as a consequence of faster transport of 
O2 from the upper stratosphere, where it is produced from CO2 photolysis, to the troposphere, 
where it is consumed by reaction with H2. We then raised and lowered the assumed stratospheric 
temperature by 20 degrees. Lowering the stratospheric temperature raised the O3 column depth 
by ~ 30%, while raising it lowered the column depth by a factor of 2. Both effects are caused by 
the increased rate of O3 destruction at high temperatures. 

We also performed one additional sensitivity test in which the planet's surface temperature 
was reduced to 240 K. This is similar to the "Snowball Earth" simulation of Liang et al. (I2006I I. 
In this case, we did not turn off volcanic outgassing of reduced gases, because as noted in the 
Introduction, this is one of the known "false positive" that could lead to arbitrarily high levels of 
atmospheric O2. Liang et al. did not predict high O2 levels, but that is because they assumed that 
photochemically produced oxidants, specifically H2O2, would accumulate in the ice. Clearly, that 
could not continue indefinitely, and so such a simulation does not represent a true steady state. In 
our simulation, with reduced volcanic gases assumed to be present, the predicted concentrations 
of O2 and O3 were both low (Table |2|. However, abiotic CH4 concentrations increased to ~ 
2.7 X 10"^ in this simulation, as a consequence of the reduced concentrations of tropospheric H2O 
and correspondingly low number densities of OH radicals. Hence, this simulation reinforces the 
thought that one should be cautious about assuming that CH4 by itself is a bioindicator, as there 
clearly are situations in which abiotic CH4 might become abundant. 

Although the sensitivity tests performed here by no means exhaust the possible sources of un- 
certainty in the photochemical model, they show that most such uncertainties lead to only modest 
changes in the results. Even a factor of 4-5 uncertainty in the predicted O2 and O3 concentra- 
tions is not signficant if the resulting column depths remain well below the level of detection. 
In general, the results are more sensitive to the boundary conditions imposed on the model, e.g., 
outgassing and escape of hydrogen, than they are to the details of the model photochemistry. 

We also generated spectra at visible and mid-infrared wavelengths for the 0.2 bar CO2 stan- 
dard early Earth case, both with and without outgassing, and we compared these to the spectrum 
of modern day Earth (Fig. |5]). In the visible-near-IR spectra (Fig. |5h), the principal differences 
between the two early Earth cases can be seen in the 0.86-0.89 fj.m, 0.99-1.03 jum, 1.15-1.20jt/m 
and 1.62-1.70 fim wavelength ranges, where the early Earth with outgassing shows more pro- 
nounced absorption from methane, as is expected given that the no-outgassing case has a ~ 10^ 
times smaller concentration of this compound (Table|2]l. For O2, although the no-outgassing case 
resulted in a 4-fold enhancement in the column depth, the total amount of O2 was still suffi- 
ciently small that there is no sign of any absorption from the O2 A-band at 0.76 jim in either 
case. Comparison with the modern Earth spectrum also shows a significant and near-identical 
reduction in O3 absorption (0.5-0.7 //m) for both early Earth cases. Again, the higher O3 column 
depth for the no-outgassing case was still too small to produce a discernible difference in the 
spectrum. The early Earth spectra also show strong CO2 absorption in a triplet of bands centered 
around 1.05 fim and a doublet of bands from 1.2-1.25 fim, features that are absent from the mod- 
ern Earth spectrum. The strongest CO2 absorption however is seen from 1.52-1.65 yum, and is 
present in the modern Earth spectrum, although at much reduced strength. 

In the mid-infrared (Fig. [Sj)), again, the differences in the two early Earth spectra are dom- 
inated by methane absorption near 7.7 yum, which is much more pronounced for the outgassing 
case. Any difference in the ozone band at 9.6 jum, however, is masked by the CO2 doubly-hot 
bands at 9.4 yum, as noted by Selsis et al. 120021 In comparison to the modern Earth spectrum, 
the early Earth MIR spectrum is dominated by the absorption of CO2 at 5.3, 7.4, 8.0, 9.4, 10.5 
and 15/vm. These bands are associated with different isotopes of CO2. The V2 fundamental of the 
dominant CO2 isotope ('^O'^C'^O) produces the strongest feature at ISfim, which has nearly 
doubled in apparent width for the 0.2 bar CO2 case, when compared to Modern Earth. The bands 
at 9.4 and 10.4 fj.m are "doubly hot bands" of the dominant isotope, transitions from upper states 
to the 2nd excited state, rather than the ground state. However, the bands near 7.4 and 8//m 
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Fig. 5. Spectra of planets with 0.2 bars of CO2 around the Sun. Standard case (red) and no out- 
gassing case (dark blue). Present Earth spectrum is plotted for comparison (dashed line). 

are not hot bands, but rather vi fundamentals for the asymmetric CO2 isotopes ('^O'^C'^O and 
'^0'^C"^0). This band is forbidden by a Fermi resonance in all symmetric CO2 isotopes, which 
includes the dominant "'O^^C'^O molecule. The flat base of the 15yum CO2 V2 absorption in the 
early Earth case is produced by the isothermal stratosphere in our model, which will also serve 
to enhance the apparent molecular absorption in the spectrum. Should the planet's actual vertical 
temperature profile deviate from this assumption, then the observed features in the MIR may be 
less pronounced, even given the same atmospheric composition. Collisionally-induced absorp- 
tion was also included in the model used to generate these spectra, but it shows relatively little 
effect on the spectrum throughout this wavelength range at these CO2 partial pressures. 
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3.2. High CO2 planetary atmospheres under high UV radiation 

In this experiment, we simulated three planetary atmospheres, each with 0.02, 0.2 and 2 bars 
of CO2, with high incoming UV radiation. The purpose was to explore the effects of high UV 
insolation on abiotically produced O2 and O3 in high-COi atmospheres. For the 0.2 bar-COi 
atmosphere under high UV, the column depths of O2 and O3 were ~3 times that of the comparable 
standard case with solar insolation (Tables |2] Fig. |6]l. These concentrations, however, are still 
very low compared with the present Earth's ozone and oxygen column depths. H2 amounts in 
the standard and high UV atmospheres were very similar (Fig.|6^), but larger photolysis rates of 
CO2, resulted in more O2 (and hence, O3) in the planet's stratosphere (see Fig.|2^). 

The ozone column depth changes very little when the CO2 rises from 0.02 bars to 0.2 bars, 
given the same amount of UV radiation received at the top of the atmosphere. But, when the CO2 
is increased to 2 bars in a high UV environment, the ozone column depth decreases to 0.76 times 
that in the planet with 0.2 bar CO2 around EK Dra. The likely reason is that CO2 photolysis 
occurred higher in the stratosphere, where densities are low and where 3-body formation of O2 
and O3 is inhibited (Fig. |2ll. In any case, the calculated changes in O2 and O3 are too small to 
produce an appreciable effect on the spectra. By contrast, the result of going from the low-UV to 
the high-UV environment, keeping the same concentration of CO2 in the planetary atmosphere, 
is to increase O2 and O3 column depths by a factor of ~3 (Table|2]l. 

Because the planet with the atmosphere containing 2-bars of CO2 had a surface T (312K) 
significantly higher than our 0.02 and 0.2 bar cases (278K), we also simulated a planet around 
EK Dra with 2 bars of CO2, but using the colder surface temperature of 278K (Fig.[T]i. This case 
was done to minimize the effect of the temperature on the ozone and oxygen destruction and to 
model an atmosphere more similar to early Mars. The O3 and O2 column depths for this case 
were less than those obtained for the hotter 2 bar-C02 planet around EK Dra, being 3.7 x 10'^ 
cm"^ and 2.7 x 10'^ cm"^, respectively. 

The spectra for these three high-UV cases (Fig.|9]to[TTJ do show differences from the standard 
model spectrum. The most obvious difference in the visible spectrum is the change in the slope of 
the spectrum between 0.5-0.7//m. This slope can be affected by both O3 absorption and Rayleigh 
scattering from atmospheric gases. Since the total O3 column depth for the 0.02, 0.2 and 2-bar 
CO2 atmospheres is nearly identical, the difference in the reflectivity in this wavelength range 
is instead caused by enhanced Rayleigh scattering, and is not the result of differences in ozone 
absorption. This scattering enhancement is due to either changes in total atmospheric pressure, 
or atmospheric composition The more massive 2-bar CO2 atmosphere has a surface pressure of 
2.9 bars, compared with 1 bar for both the 0.02 and 0.2 bar CO2 cases, and its large enhancement 
in the Rayleigh scattering slope between 0.5-0.7 jjm is due primarily to its larger atmospheric 
pressure. However, the 0.02 and 0.2 bar CO2 atmosphere cases also show a small difference in 
their Rayleigh scattering slopes, even though they have identical atmospheric pressures. This is 
due instead to the difference in the composition of their atmospheres, as CO2 is a more efficient 
Rayleigh scatterer than nitrogen, and enhances the Rayleigh scattering slope for the 0.2 bar CO2 
case over the 0.02 bar case. In addition to the Rayleigh scattering slope, CO2 absorption fea- 
tures are also seen strongly in the spectrum and reveal spectral sensitivity to atmospheric CO2. 
The triplet of absorption bands near l.OSyum produces no appreciable absorption for the 0.02 bar 
case, but can be distiguished as a triplet at high resolution (R ~ 140 and above) for both the 0.2 
and 2-bar CO2 cases (Fig. |9] - inset). At resolutions of R ~50 and greater, the triplet structure 
cannot be clearly discerned, but the ability to distinguish between the 2-bar CO2 and the other 
two CO2 concentration cases at 1.05yum is maintained (Fig. [TOb. assuming high S/N measure- 
ments. The CO2 bands at 1.206 and 1.221/vm become progressively stronger as the atmospheric 
CO2 abundance is increased and can also be distinguished as two features at R ~ 140, but still 
retain sensitivity to atmospheric CO2 abundance at resolutions as low as R ~ 50. A similar trend 
is shown for the stronger CO2 bands near 1 .6 jjm, however in this case, the 2-bar CO2 atmo- 
sphere is so heavily absorbing that the planetary spectrum becomes black throughout much of 
the 1.52-1.7yum region, with the exception of a relatively transparent, narrow atmospheric win- 
dow at 1.55//m. 

Finally, marked differences in the strengths of the water bands near 0.94/im, 1.14yum and 
1.4yum are observed for the 2-bar CO2 case. This enhanced absorption is the result of increased 
H2O abundances for the 2-bar CO2 case, as compared to the two lower-C02 cases (FiglT]i, due 
to the higher surface temperature. The model tropospheres assume a fixed distribution of relative 



Segura et al.: Abiotic O2 and O3 in higti CO2 atmospheres 



15 



a) 





60 r 




50 \ 


E 


40 f 








30 \ 


< 


20 \ 






llllllll 



Present solar LTV // 

// 



Higher UV 



10"^' 10-'^ 



10 10 10 
O^and mixing ratios 



10" 



60 r 
50 \ 
40 \ 
30 f 
20 I- 
10 I- 

oL 



10^ 



Present solar UV 



I \ 
} 

/ 

/ 



Higher UVj 



10' 



10^ 



10^ 



10^ 



Ozone number density (cm ) 



Fig. 6. O2 mixing ratios (a) and O3 number densities (b) for atmospheres with 0.2 bars of CO2 
under present UV solar radiation and high UV flux from a solar-like young star 



humidity; hence, H2O abundances increase exponentially with increasing surface temperature, 
following the corresponding increase in the saturation vapor pressure. The water absorption fea- 
tures can be identified, and different concentrations discriminated between, at resolutions as low 
as R ~ 20. However the ability to separate other species from the water vapor (such as the 1 .2yum 
CO2 absorption feature) at R ~ 20 becomes extremely difficult, instead requiring resolutions 
closer to R ~ 50 - 70 to do so, even for high S/N. 

In the MIR, the spectra are dominated again by absorption from CO2, and these early Earth 
spectra are not at all similar to that of modern Earth (Fig.fTTh). In particular, for the 2-bar CO2 
case, the spectrum is heavily distorted by the CO2 absorption, which narrows the Earth's normal 
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atmospheric window region between 8-13yum down to only 8-9jum. The effects of collisionally- 
induced absorption are more pronounced here for the 2-bar CO2 case, and manifest themselves 
in this wavelength range primarily as a broadening of the CO2 isotope bands, enhancing the wing 
absorption near 7.4 and 8.0 ^im. While there is some theoretical justification for CO2 pressure- 
induced continuum absorption within the 10 [im and 15 [im bands, existing measurements show 
that this absorption is apparently very weak compared to the allowed transitions at these wave- 
lengths. We are not aware of any quantitative estimates of its intensity. For example, recent 
high-pressure (200 atmosphere) measurements by Hartmann's group (c.f. Niro et al . 120041 1 show 
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Fig. 8. O2 mixing ratios (a) and O3 number densities (b) for atmospheres with different amounts 
of CO2 under high UV flux from a solar-Hke young star 



~ 10% residuals with respect to their best line mixing model in simulations of absorption at pres- 
sures at 200 atmospheres. Some of these residuals might be associated with collision-induced ab- 
sorption, which is neglected in their model. In any case, if this absorption scales as p^, it should 
be lO'* times weaker for a 2-bar atmosphere, and would therefore have a negligible effect on the 
modeled spectrum within the context of the work reported here. 

In all three cases, methane absorption at 7.7jum is greatly enhanced over that of modern Earth, 
given that higher concentrations of this compound are present in these I0W-O2 cases (Table |2]i. 
Looking at the MIR brightness temperatures for these three cases (FigfTTb). the atmospheric 
window in the 2-bar CO2 case exhibits an enhanced brightness temperature of 287K when com- 
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Fig. 9. Visible-NIR spectra of planets with 0.02 (dark blue), 0.2 (red) and 2 bars (black) of CO2 
around EK Dra. 

pared with the 272K brightness temperature for the other two cases. This is due to the fact that 
the atmosphere is considerably hotter overall, but note that the surface temperature for the 2-bar 
CO2 case is 317K, and for the other two, 278K. In all cases, even though no clouds were present, 
it was not possible to retrieve the surface temperature of the planet by measuring the brightness 
temperature in the atmospheric window. This is due to the fact that the atmospheric window is 
not completely transparent, but includes weak absorption from water vapor in an Earth-like atmo- 
sphere, and potentially other absorbers in a non-Earth-like atmosphere. The water vapor, which 
is much enhanced in the 2-bar CO2 case has effectively truncated the atmospheric column at a 
level higher and cooler than the surface. For the 2-bar atmosphere the 30K discrepancy between 
the actual surface and the measured brightness temperatures corresponds to a sampled altitude of 
~ 4km above the surface. 

Direct comparison between spectra generated for the 0.2 bar CO2 planetary atmospheres for 
high-UV (EK Dra) and low-UV (Sun) cases show very few differences, except for changes in 
the strength of the methane absorption features for both the vis-NIR and MIR. This is due to 
the UV-induced removal of CH4 from the EK Dra planet, which reduces the overall CH4 mixing 
ratio by a factor of ~3 for the high-UV case, as can be seen in Table |2] 

Finally, the modeled strength of the O2 A-band feature in the visible, for all the planetary 
cases considered here (FigfTSli demonstrates that the generation of a detectable false positive for 
oxygen is highly unlikely, even with very high spectroscopic resolution and S/N. 

4. Discussion 

Our results show that for a habitable planet, rainout of oxidized species onto a reduced surface 
will inhibit the generation of an abiotic "false positive" signal with strong O2 or O3 absorption. 
As summarized in Figure [12] this result is robust to both high stellar UV output and high at- 
mospheric CO2 abundances, as abiotic O2 and O3 are not visible for any of our cases, even in 
the extreme case of a planet with an atmosphere containing 2-bars of CO2 and in orbit around a 
high-UV star. We also explored the sensitivity of atmospheric constituents and their detectability 
to volcanic outgassing rates, incident UV flux, and atmospheric composition. We found that the 
most significant spectral effects of either turning off volcanic outgassing, or increasing UV flux, 
was to noticeably decrease the observed methane absorption. Increased stellar UV flux had little 
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Fig. 10. Spectra for the 0.02, 0.2 and 2-bar CO2 cases convolved with a triangular slit function 
and degraded to a constant wavenumber resolution equivalent to a)R ~140, and b) R ~50. 



observable effect on the rest of the planetary spectrum, even though modest enhancements were 
seen in the atmospheric O3 and O2 column depths. 

Our results contradict previous models (Selsis et al. I2002I I that predicted high abiotic abun- 
dances of O2 and O3 in planets with COi-rich atmospheres. As explained earlier, these models 
failed to balance the H2 budget because they ignored the rainout of soluble oxidized and reduced 
compounds. Including these terms-particularly the rainout of oxidants such as H2O2, H2SO4, 
and HNO, ensures that an abiotic atmosphere on an Earth-like planet should contain substantial 
concentrations of H2 and only minimal amounts of O2 and O3. 
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Fig. 11. a) MIR spectra for present Earth (green), our standard case (blue, 0.2 CO2 bar, present 
solar UV) and 2 bars (black) of CO2 around EK Dra. b) MIR spectra of planets with 0.02 (dark 
blue), 0.2 (red) and 2 bars (black) of CO2 around EK Dra. 



That said, there still may be planets within the habitable zones of their parent stars that do 
not have active hydrological cycle for which calculations similar to those of Selsis et al. may be 
valid and on which abiotic O2 and O3 could conceivably build up. For example, if a planet in 
the habitable zone lacked significant liquid water, perhaps because it formed with no appreciable 
water inventory, as has been modeled for teiTestrials in planetary systems with eccentric Jovians 
(Raymond et al. 120041 1. then our entire analysis of the hydrogen budget would be invalid. H2 
would not be produced by reaction of photochemically produced oxidants with reduced species 
in the crust and ocean, nor would it be likely to be outgassed from volcanoes. This would be 
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Fig, 12. O2 A-band for all the simulated atmospheres. 



compensated by a complete lack of hydrogen escape to space, thereby eliminating the net source 
of O2 identified earlier in this paper. Some abiotic O2 would presumably accumulate in the atmo- 
sphere due to photolysis of CO2, although we have not attempted to quantify O2 production for 
dessicated planets in this paper. However, abiotic production of O2 on such a water-free planet 
could be relatively easy identified by the complete lack of H2O bands in the planet's visible and 
MIR spectrum, which would indicate that it was unlikely to be an abode for Earth-like, water- 
and-carbon-based life. 

Another exception to our analysis would be a frozen, "Snowball Earth" type planet within the 
habitable zone of its star. Earth itself appears to have undergone global glaciation at least 3 times 
during its history, most recently in the Neoproterozoic, at ~ 610 Myr and 720 Myr ago (Hoffman 
et al. I1998I I. On such a frozen planet, as on the dry planet described above, photochemically 
produced oxidants would not react with reduced materials in the planet's crust and ocean. The 
planet would thus be similar to the Mars-like planet mentioned in the introduction for which 
atmospheric O2 could, in theory, continue to accumulate leading to a possible "false positive" for 
life. 

Two things mitigate against this posing a significant problem for life detection by TPF or 
Darwin. One is that Snowball Earth episodes on Earth are thought to be relatively short-lived, lO's 
of millions of years or less, because they are eventually terminated by buildup of volcanic CO2 in 
the atmosphere (Hoffman et al. 1998 ). 10 Myr represents only about 0.2% of the Earth's geologic 
history, and the odds of finding an Earth-like planet in such a state in our Solar neighborhood are 
relatively small. Second, it should be possible to identify whether a planet is or not in a Snowball 
Earth-type phase, either by looking for the presence of water ice or water vapor bands in their 
spectrum, or by follow-up measurements of polarization to characterize the planetary surface and 
search for the presence of liquid water. 

To quantify the effect of decreasing global surface temperature on the detectability of water 
vapor we ran a series of models similar to our standard, 0.2 bar CO2 case, but for successively 
lower surface temperatures. The corresponding synthetic spectra generated from these environ- 
ments are shown in Figs. [T3] and fT4l Model calculations of Snowball Earth (Hyde et al. I2000I I 
show that the surface temperature drops precipitously from ~ 275 K to 230 K when global glacia- 
tion sets in. For an atmosphere with Earth's relative humidity distribution, Fig.[T3]shows that the 
water vapor bands in the visible and near-IR become extremely weak as this transition occurs 
with a 238K surface temperature planet showing only ~ 5% of the absorption strength in the 
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Fig. 13. Albedo of planets with 0.2 bars of CO2 around the Sun with different surface tempera- 
tures. 



1.13;um water band when compared to the feature seen for the 278K surface temperature planet. 
However, with high S/N and a spectral resolution greater than R -50, the 1.13/im feature could 
be detected at temperatures greater than 238K. It must be stressed that S/N and resolution are 
trade offs, as lower resolutions can require higher S/N to detect a given feature. In the MIR (Fig. 
[T4l i the changes in the spectrum with decreasing surface temperature are seen as an interplay 
between lower water vapor amounts and the atmospheric temperature. The declining planetary 
atmospheric temperatures with surface temperature are seen in these cloud-free cases in the at- 
mospheric window region between 8-13/im. Although note that even in these cloud-free cases, 
the peak temperature detected in the atmospheric window at 9fim is not the surface temperature, 
being several degrees cooler due to the presence of water vapor, which results in atmospheric 
temperatures sampled at levels higher than the surface (Fig.[T4b). If clouds were present, the in- 
ferred temperature would typically be even colder. The declining water vapor amounts are seen in 
the difference in the brightness temperature between the 5.5-7.0/im water vapor region and the 8- 
13/im window region. Correct interpretation of these spectra, however, would require generation 
of synthetic spectra that could be used to disentangle the combined effects of water and temper- 
ature on the spectrum. Nonetheless, as was the case more straightforwardly for the visible-NIR 
spectrum, it may be possible to detect the significantly lower water vapor amount, given instru- 
mentation that was sensitive to the 6.3fim water vapor band. In this case, a resolution greater than 
R ~ 20 in the MIR might be sufficient, with high S/N. However, we add the caveat that the pres- 
ence of clouds, either planet-wide, or covering a smaller fraction of the visible planetary disk, 
could skew the interpretation to colder temperatures and lower water abundances than are present 
at the planet's surface. For these cases, searching for abundant water-soluble atmospheric gases, 
such as SO2 might serve as a secondary indicator that we are observing a planet without a sur- 
face ocean, rather than one with clouds obscuration. Multiple samples of the planet's atmosphere 
over time may also help to determine whether clouds are affecting the observed temperatures and 
water abundances. 

Another way to test for the presence of liquid water on an exoplanet's surface is to measure 
the polarization of light reflected from its surface. The specular "glint" produced by the smooth 
ocean surface should be strongly polarized, at least at some phase angles (McCullough 120071 
Stam et al. i2006. West et al. 12005 1 Williams & Gaidos 2004 ). Separation of different polarizations 
is a feature of some suggested techniques for coronagraphic starlight suppression. Even if this 
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Fig. 14. MIR spectra of planets with 0.2 bars of CO2 around the Sun with different surface tem- 
peratures.a) Radiance b) Brightness temperature 



measurement was not performed by an initial TPF-C mission, it could presumably be done by 
some follow-up mission as a way of verifying the presence of surface liquid water (or at least 
of some surface liquid). A positive detection of either O2 or O3 in an exoplanet's atmosphere 
would certainly warrant further investigation to be sure that the liquid water required for life was 
present at the planet's surface. 
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4.1. Applicability to other planetary systems 

For these simulations, we used EK Dra, a young, UV-active star with a calculated age of only 
10^ years. If we were considering only whether this star was likely to harbor habitable planets, 
then this would be a poor choice, as planets around a star of this age would still be undergoing 
large, ocean-vaporizing impacts (Sleep et al.'1989). We would not anticipate seeing planet-wide 
biosignatures at this early and hostile stage, and so would also be unlikely to mistake a false 
positive as a biosignature. Nonetheless, EK Dra is an excellent example for our simulations, 
as it most realistically represents the approximate upper limit on the UV radiation likely to be 
encountered by a planet around a Sun-like star. 

When the first studies of the UV environment for early Earth were published, the only known 
examples of young Sun-like stars were the T Tauri type stars. lUE observations showed that T 
Tauri UV fluxes were as much as 10^ times higher than the present solar UV radiation. Based 
on those observations and stellar evolution models, those early studies concluded that the UV 
enhancement for a 10^ year old. Sun-like star, should be of the order of 10 times the present 
solar UV flux (Canuto et al. [19821 Zanhle & Walker '1982'). However, we believe that this is an 
overestimate. EK Dra, at 10* years, and one of the most UV active stars known (Ribas et al l2005l l. 
does exhibit fluxes 10 times that of the Sun at wavelengths shortward of 160nm. However, when 
this flux is integrated over the wavelength range most relevant to atmospheric chemistry (A < 200 
nm), it is only ~ 4.5 times larger than that of the Sun, on average. More importantly for the 
atmospheric chemistry, the continuum photon flux (rather than the energy flux) at A < 200 nm 
is just ~ 3.6 times higher in EK Dra when compared to the solar spectrum (c.f. [3]) (Note that 
our continuum calculations do not include the Lya line which we calculated to be ~ 10 times 
larger than for the present Sun, and which was included in all our modeling simulations). These 
considerations suggest that our results can be interpreted as an upper limit on detectability for 
high-C02 atmospheres around young, UV-active stars that are either Sun-like, or cooler than 
our Sun. If we consider stars with higher-UV flux than our Sun, e.g. the F stars, it is also highly 
likely that our results are applicable. An F2V star has an integrated UV flux that is approximately 
2.5 times that of EK Dra (Segura et al. l2003l l. An examination of our Table 2 indicated that the 
increase in O2 and O3 column depth in a high-COi atmosphere as a result of EK Dra's 4.5 times 
higher flux is relatively modest, being 3.6 and 2.9 times respectively, and this is for column 
depth values that are still 10^ times lower than those likely to be detectable for a planet around an 
F2V star (see Segura et al. 120031 Table 1 and Fig 14). It is therefore unlikely that the additional 
increase of 2.5 in the UV flux from an F star would result in detectable O2 and O3 for the high- 
CO2 planets modeled here. We therefore conclude that the results presented here are an upper 
limit for the majority of parent stars likely to be considered as TPF/Darwin targets. 

4.2. Effect of clouds on our conclusions 

All the simulations shown here were run for clear-sky atmospheres, with no appreciable hazes 
or cloud layers, and isothermal stratospheres. These cases shown here are therefore the most 
sensitive for detecting O2 and O3 in our planetary spectrum, as they allow a longer pathlength 
through the atmosphere that is not truncated by clouds. Additionally, the contrast between the 
hot surface and lower atmosphere, and the cool, isothermal stratosphere provides for maximum 
detectability in the mid-infrared. However, given that both the O3 and O2 concentrations peak 
above the tropopause in all our model cases, and have significantly smaller abundances below 
that, a cloud at the tropopause would not significantly change the net column abundance of O2 
and O3, and hence their detectability. In the case of the CO2 however, which is evenly mixed 
through our atmospheres, a cloud at the tropopause would significantly truncate the observed 
column, and would result in much reduced absorption depths for CO2 in the visible, and would 
also dramatically reduce the CO2 absorption in the mid-infrared, as the cooler clouds tops would 
obscure the hotter surface and lower atmosphere and provide a reduced temperature diff'erence 
with the stratosphere. The brightness temperatures of the observed emitting layer would also 
be reduced to that of the cloud top temperatures, and would preclude our ability to determine 
the planet's surface temperature. This case can be demonstrated by comparing the brightness 
temperatures for the 2-bar CO2 atmosphere and a modeled spectrum of the planet Venus (Fig.fTsTl. 
Even though Venus is known to have a 93bar CO2 atmosphere and a 730K surface temperature. 
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Fig. 15. Spectra of the 2 bars of CO2 planet around EK Dra and Venus. 



its CO2 absorption is significantly weaker than seen in our cases, and its maximum brightness 
temperature in the MIR is only 250 K, due to its planet-wide cloud deck. 

4.3. Implications for extrasolar planet characterization missions 

TPF-C, TPF-I, and Darwin will all have the opportunity to study a wide variety of planetary 
systems. If terrestrial planets are found around these stars, it is highly likely that we will see 
the atmospheres of planets, like our own Solar System examples of Venus, Mars and the early 
Earth, that are dominated by CO2. Therefore, understanding the appearance of C02-dominated 
atmospheres, and their abiUty to generate false positives is an important step in characterizing 
extrasolar terrestrial planets that are not strictly Earth-like. In addition to disproving the likeli- 
hood of significant abiotic sources of O2 and O3 in a C02-dominated atmosphere, we have also 
demonstrated that C02-dominated habitable atmospheres exhibit a wealth of spectral features. 
The multiple CO2 absorption bands throughout the visible-NIR range accessible to an optical 
coronograph (0.5-1.7/im) and an MIR interferometer (5-25yum) provide differing sensitivity to a 
large range of CO2 abundances. The most sensitive band in the visible-NIR is near 1 .62/im, and 
produces a potentially detectable feature, even for the Earth's present 330ppm abundance. For 
atmospheres with higher CO2 abundances (0.2-2bars of CO2), the less sensitive 1.05;um bands, 
which are weak in the Earth's atmosphere, provide strong absorption features, and are better fea- 
tures for potential quantification of the planet's atmospheric CO2 provided the spectral resolution 
is above 50. In the MIR, the CO2 15;um feature provides the strongest absorption for CO2 and 
is the most detectable feature in a terrestrial planet atmospheres. In addition, we have shown 
that in the MIR the isotopes of CO2, produce strong features. This is arguably one of the most 
detectable isotope features in a terrestrial planet spectrum, even at relatively poor spectral resolu- 
tion (R ~ 20) and signal to noise. These features may allow us to get a first measure of O isotopic 
composition of an early-Earth-Uke or C02-dominated terrestrial planet. As Ughter isotopes tend 
to evaporate more readily, and heavier isotopes are rained out more efficiently, the atmospheric 
O isotope inventory may in turn provide secondary clues to the presence and degree of activity 
in a planet's hydrological cycle. 
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5. Conclusions 

Our results show that for a habitable planet with abundant Uquid water, rainout of oxidized 
species onto a reduced surface will keep atmospheric H2 sufficiently high, that the generation 
of an abiotic "false positive" signal with strong O2 or O3 absorption is unlikely, and does not 
pose a significant hazard for habitable planets observed with TPF or Darwin. This result is ro- 
bust to realistic Umits on the UV output of the star, and the abundance of CO2 in the planet's 
atmosphere. Dense CO2 atmospheres of rocky planets may also display a range of distinctive 
spectral features in both the visible and MIR that allow us to at least qualitatively determine the 
prevalence of CO2 in the atmosphere, and potentially measure atmospheric isotopic ratios for 
oxygen. 
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